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Abstract 

Context. Efficient particle acceleration can modify the structure of supernova remnants. In this context we present the results of the 
combined analysis of the XMM-Newton EPIC archive observations of SN 1006. 

Aims. We aim at describing the spatial distribution of the physical and chemical properties of the X-ray emitting plasma at the shock 
front. We investigate the contribution of thermal and non-thermal emission to the X-ray spectrum at the rim of the remnant to study 
how the acceleration processes affect the X-ray emitting plasma. 

Methods. We perform a spatially resolved spectral analysis on a set of regions covering the whole rim of the shell and we exploit the 
results of the spectral analysis to produce a count-rate image of the "pure" thermal emission of SN 1006 in the 0.5-0.8 keV energy 
band (subtracting the non-thermal contribution). This image significantly differs from the total image in the same band, specially near 
the bright limbs. 

Results. We find that thermal X-ray emission can be associated with the ejecta and we study the azimuthal variation of their physical 
and chemical properties finding anisotropies in temperature and chemical composition. Thanks to our thermal image we trace the 
position of the contact discontinuity over the whole shell and we compare it with that expected from 3-D MHD models of SNRs with 
unmodified shock. 

Conclusions. We conclude that the shock is modified everywhere in the rim and that the aspect angle between the interstellar magnetic 
field and the line of sight is significantly lower than 90° . 

Key words. X-rays: ISM - ISM: supernova remnants - ISM: individual object: SN 1006 



^ 1. Introduction (iBerezhko & Ellisonll999llDecourchelle et alj2000llBlasil2002l 

5—1 and references therein). All these effects modify the spectral 

The diffusive Fermi acceleration process in Supernova Remnant properties of the thermal emission. Therefore, both thermal and 

(SNR) shocks can account for the observed spectrum of galactic no^.thermal emission can provide important information on the 

cosmic rays at least up to the knee at 3 x 10 ^ eV and even be- acceleration processes and can be useful to verify the expected 

yond (. Berezhko&Volk2007, Blandford & Eichleii a987). The ^f^^^^^ ^^^^^ modification and to study the open questions 

AS CA detection of synchrotron X-ray emission in SN 1006 has ^^^jj ^^^^^ j^^ate, for example: what are the factors that influ- 

shown that electrons can be accelerated up to very high energy j^e fraction of particles injected in the acceleration pro- 

at the shock front, thus providing a direct proof for the connec- ^^^^^ ^j^^^ ^^eir maximum energy? What is the relationship 

tionbetween SNRs and cosmic ray acceleration d Koyama et al. | ^^^^^^^ the upstream/downstream magnetic fields and the ac- 

| 1^^3|) . celeration processes? 

The loss of energy "deposited" in the production of rela- 
tivistic particles and their non-Unear back-reaction on the back- To this end SN 1006 is a privileged target. Its rather simple 
ground plasma are predicted to strongly affect the structure of morphology, the exact knowledge of its age, and its evolution 
the remnant, by increasing the shock compression ratio well in a low density uniform medium (because of its high galactic 
over the Rankine-Hugoniot limit and decreasing both the post- latitude) makes it an ideal laboratory for these studies. The bilat- 
shock temperature and the geometric extension of the region eral morphology of SN 1006 shows two opposed radio and X-ray 
between the main shock front and the contact discontinuity bright limbs characterized by non-thermal emission separated by 
a region with low surface brightness. The synchrotron limbs re- 
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face brightness (see, for exampl e . iRothenflug et al.ll2004l) . As 
pointed out by IRothenflug et alj (|2004|) . the very low surface 
brightness above 2 keV in the interior strongly suggests that the 
ambient magnetic field is parallel to the shock velocity in the 
bright limbs (polar-cap scenario), in the framework of the quasi- 
parallel model for the injection (i. e. the injection efficiency of 
the accelerated particles increases with decreasing the obliquity 
angle, t/r, betwee n the magnetic field and the shock velocity). 
IRothenflug et alJ (2004) studied the azimuthal variations of the 
spectral properties of the non-thermal emission over the whole 
shell and its radial profile from the center to the north-eastern 
bright limb. They found that the roll-off frequency of the syn- 
chrotron emission is much larger in the bright non-thermal limbs 
than in the center and elsewhere in the rim. These results are 
difficult to explain in the polar cap scenario, but are quite un- 
expected also in the "equatorial belt" scenario (i. e. the ambient 
magnetic field. Bo, is oriented from South-East to North- West 
and the injection efficiency is larger when ifi ~ 90°), which pre- 
dicts variations of a factor ~ 2, that is mu c h less than the ob- 
served factor > 10. Recently, ' Petruk et al.l (l2008l) have shown 
that, by assuming a uniform ambient magnetic field in the polar- 
cap scenario, SN 1006 should be a centrally filled SNR in radio 
to reproduce the observed azimuthal radio profiles and this is 
contrary to what is observed. The determination of the geometry 
of the magnetic field in SN 1006 and its effects on the accel- 
eration process are therefore an open problem and the study of 
thermal emission may provide useful indications, since it can be 
strongly affected by shock acceleration processes. 

The X-ray thermal emission is, in comparison to the non- 
thermal, more un iformly distributed over the whole remnant 
(IRothenflug et al.ll200 4). Despite the importance of thermal X- 
rays for the study of shock acceleration processes, the physical 
origin of the thermal emission is still uncertain and the spatial 
distribution of its properties has not been studied yet. 

Several studies have shown that thermal X-rays in SN 1006 
cannot be u niquely associate d with solar-abundance shocked 
plasma (e. g. iDver et al.ll2004l) . X-ray emitting ejecta have been 
obser ved in th e north-western and south-eastern rims of the rem- 
nant (lAcero e t al. 2007) and even Fe K lines emission has been 
recently detected in the inter ior of the shell, where iron over- 
abundance has been revealed (lYamaguchi et al.l2008! ). The pres- 
ence of significant X-ray emission from the shocked interstellar 
medium (ISM) and the value of its temperature are indeed still 
controversial. Acero et al. (2007) modelled the X-ray emission 
of their spectral regions with one non-thermal component plus 
two thermal components (ejecta for the soft emission and ISM 
for the hard emission). However, they found that the shocked 
ISM component (with kT/sM ~ 1.5-2 keV) is almost not needed 
from a statistical point of view, since the quality of the fits does 
not ch ange by associating the thermal emission only with the 
electa. lYamaguchi et al.l(l2008l) . instead, used three thermal com- 
ponents to model the SUZAKU thermal emission of very large 
regions of SN 1006, by associating the soft component with the 
shocked ISM {kTjsM ~ 0.5 keV), and the hot components with 
the ejecta. Nevertheless they cannot exclude that the O line com- 
plexes, dominating the soft component, originate in the ejecta. 

A first approach to study the shock modification ef- 
fects through the X-ray th ermal emission was pursued by 
ICassam-Chenai et al.l (l2008l) (hereafter COS). They focussed on 
the south-eastern quadrant of the shell and they assumed the 
emission in the 0.5 - 0.8 keV band to be completely associ- 
ated with the ejecta. They measured the ratio, rgw/CD, of radii 
between the blast wave shock front (BW) and the contact dis- 
continuity (CD) in several points, finding very low values with 



a minimum rgw/CD ~ 1 near the non-thermal limbs. The az- 
imuthal trend of rgw/CD seems to be in agreement with the polar 
cap scenario, although the values of rgw/CD are too low and not 
consistent with that model. 

Since a spatially resolved study of the thermal emission at 
the rim of SN 1006 has not been performed yet, here we present 
the analysis of archive XMM-Newton observations of SN 1006. 
We focus on the rim of the remnant to study the azimuthal vari- 
ations of thermal and non-thermal X-ray emission immediately 
behind the main shock front. We aim at obtaining detailed infor- 
mation on the nature of the thermal emission and to study its link 
with the shock acceleration processes. We also use a Chandra 
archive observation of the north-eastern limb of SN 1006 taking 
advantage of its larger spatial resolution to test our conclusions. 
The paper is organized as follows: the data and the data anal- 
ysis procedure are described in Sect. |2] in Sect. [3] we present 
the X-ray results in terms of spatially resolved spectral analysis 
and image analysis. Our results are discussed in Sect. |4| and our 
conclusions are summarized in Sect.|5] 

2. Data processing 

We consider all the XMM-Newton EPIC archive observa- 
tions of SN 1006 available, but we do not present here the 
observation 0111090601 since it is strongly affected by pro- 
ton flares and its cleaned exposure time is quite low 1 ks). 
All the observations that we analyze were performed with the 
Med ium filter by using the Full Frame Mode for the MOS cam- 
eras dTurner e t al.' 2001 ) and the Ext ended Full Frame Mode for 
the pn camera ( St riider et al.l l200lb . The relevant information 
about the XMM-Newton observations presented here are sum- 
marized in Table [1] We process the data by using the Science 
Analysis System (SAS V7.1). Light curves, images, and spec- 
tra, were created by selecting events with PATTERN<12 for the 
MOS cameras, PATTERN<4 for the pn camera, and FLAG=0 
for both. To reduce the contamination by soft proton flares, we 
have screened the original event files using the sigma-clipping 
algorithm suggested by Snowden & Kuntz (2007). The screened 
exposure times are given in Table [Tj 

We produce the images of the whole remnant by superpos- 
ing (we use the E MOSAIC task) the MOSl, MOS2, and pn 
images of the pointings shown in Table [Tj We have taken into 
account the differences between MOS and pn effective areas 
and the mosaiced images are MOS -equivalent. All the images 
are background-subtracted, vignetting-corrected, and adaptively 
smoothed (with the task AS MOOT H). The exposure and vi- 
gnetting corrections were performed by dividing the count im- 
ages by the corresponding superposed exposure maps (obtained 
with the task EEXPMAP). 

Spectral analysis was performed in the energy band 0.5 - 5 
keV using XSPEC. To check the validity and the robustness of 
our result we decide to adopt two different procedures for the 
spectral extraction and for the background subtraction. 
1. We use the event files generated with SAS. We generate the 
Ancillary Response Files with the SAS ARFGEN task, and we 
use the EVIGWEIGHT task ( Arnaud et al. 2001) to correct vi- 
gnetting effects. The contribution of the background is derived 
from the high statistics background event files discussed by 
ICarter & ReadI (l2007l) . We subtract the background contribution 
from the same region positions on the CCD (i. e. in the detector 
coordinates), in order to take into account the inhomogeneous 
response of the instruments across the field of view. For each of 
the spectral regions presented in Sect. [3] we extract the spectrum 
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Table 1. Relevant information about the data. 



OBSJD 



Date 



MOS-pn f„p (ks)* 



RA, J2000 
15"03'"50' 
15''03"'30' 
15''01'"51' 
15''04'"07' 
15''03"'30' 
15''02'"35' 
15''03"'35' 



DEC, J2000 



0111090101 
0111090601 
0077340101 
0077340201 
0143980201 
0202590101 
0306660101 



2000- 08-20 

2001- 08-08 
2001-08-10 
2001-08-10 

2003- 08-14 

2004- 02-10 

2005- 08-21 



7-3 
7-3 
30-21 
24-19 
16-11 
26-17 
11-4 



-4r47'00" 
-42°01'00" 
_4r49'00" 
-41°52'32" 
_41°48'12" 
-42°04'37" 
-42° 04' 20" 



■ Unscreened/Screened exposure time. 



from the corresponding event file with the largest exposure time. 
Spectra were rebinned to achieve a signal-to-noise ratio per bin 
> 5cr and the fittings were performed simultaneously on both 
MOS spectra and on the pn spectrum. The reported errors are at 
90% confidence. 

2. We use the Extended Sources Analysis Software (XMM- 
ESAS V2.0) which allows us to subtract instrumental and cos- 
mic backgrounds separately. XMM-ESAS models instrumental 
background from the "first principles", using filter-wheel closed 
data and data from the unexposed corners of archived observa- 
tions. The cosmic background is then modeled in XSPEC after 
the instrumental background subtraction. The standard XMM- 
ESAS method of background subtraction has some difficulties, 
namely the large number of parameters and lack of suitable 
model for the Al K„ (~ 1.49 keV) and Si K„ (~ 1.75 keV) in- 
strumental lines. To overcome these difficulties, we develop a 
procedure which is similar to the "do uble background subtrac- 
tion" technique ( Read & Ponmanll20b3i) . In this procedure, for 
each observation we produce spectra and response files from the 
"background" region, which lies out of SN 1006. Then we model 
the "background" spectra in XSPEC, and keep these "back- 
ground" parameters fixed when fitting "source" spectra after- 
words To improve our knowledge of local cosmic background 
parameters, we add to them the RASS background datcQ as sug- 
gested in XMM-ESAS manual. The obtained best-fit parameters, 
describing cosmic background, remaining soft proton contribu- 
tion, and instrumental lines, are then fixed and used to model 
the parameters of "source" regions, after proper rescaling. To 
prepare the EPIC MOS event lists we use the XMM-ESAS 
script MOS-FILTER. We obtain the MOSl and MOS2 spectra 
and construct the corresponding background using the XMM- 
ESAS scripts MOS-SPECTRA and XMM-BACK, respectively. 
Finally, we group the spectra and the corresponding background 
to have at least 25 counts per bin. 

We here present the results obtained with method 1, since 
it allows us to use also the pn data. We verified that the results 
obtained are perfectly consistent with that produced with method 
2. 

We also use a Chandra archive observation of SN 1006 to 
perform the test described in Sect. 13.21 In particular, we use the 
observation of the NE limb of SN1006 with ID 732, taken on 7 
Jul 2000, which has an exposure time of 69 ks. 

3. The data analysis 

3.1. Spatially resolved spectral analysis 

Figure [T] shows the mosaiced EPIC images of SN 1006 in the 
0.8 - 2 keV band, 0.5 - 0.8 keV band ("oxygen band", central 



' They were obtained using HEASARC background tool, 
|http : //heasarc . nasa . gov/cgi-bin/Tools/xraybg/xraybg . pi , 



panel), and in the 2 - 4.5 keV band ("hard band", lower panel). 
Since we aim at studying the effects of the acceleration process, 
we select our spectral regions at the border of the shell. The set 
of 30 regions selected for the spatially resolved spectral analysis 
is shown (in white) in the upper panel of Fig. [T] All the regions 
cover the same area in the plane of the sky and they extend in 
the radial direction for 1.5' (i. e. ~ 10% of the SNR radius), 
corresponding to ~ 1 pc at 2.2 kpc. Our approach is similar to 
that followed by iRothenflug et al.l (l2004l) . although we here fo- 
cus also on the thermal emission. 

There are striking diff'erences among the 30 spectra. Figure 
12] shows the pn and MOS spectra extracted from region 29 
(upper panel), where several emission line complexes are vis- 
ible, and, for comparison, from region 23 (lower panel), lo- 
cated in the bright north-eastern non-thermal limb, where the 
emission is completely featureless. Despite these differences, we 
plan to describe all the spectra with a unified model, so as to 
explain the diff'erent spectral properties in terms of azimuthal 
variations of the best-fit parameters. We then use an isother- 
mal model of optically thin plasma in non-equilibrium of ion- 
ization with a linear distribution of ionization timescale ver- 
sus emission mea sure and with fre e abundances (VPSHOCK 
model in XSPEC. iBorkowski et"aLl l2b01) to describe the ther- 
mal component, plus a synchrotron emission from an electron 
power law with exponential cut-ofF (SRCUT model in XPSEC, 
[Reynolds & Keohane|[T999l) to model the non-thermal compo- 
nent. We verify that the quality of the fit does not improve 
by adding another thermal component. Moreover, an additional 
thermal component introduce too many free parameters (consid- 
ering the available statistics) and this generates entanglements 
between the best-fit values, thus determining large errors and 
useless results. We therefore consider only one thermal com- 
ponent. Our model of the thermal emission is ther efore diff'er- 
ent fro m the two-temperatures model adopted by lAcero et al.l 
(I2007h for the north-eastern and south-western regions of the 
shell an d it is simpler tha n the three-temperatures model used 
by Ya maguchi et al.l ( l2008l) . where much larger and less uniform 
extraction regions are considered. 

We constrain the normalization of the non-thermal compo- 
nent (i. e. its flux at 1 GHz) by extracting ffie radio flux from 
the same regions defined in the X-ray map. To this end we use 
a radio image obtained from interferometric VLA observations 
with the addition of single dish data. This procedure allowed 
us to produce the highest angular resolution radio image of SN 
1006, together wiffi an accurate recovery of the flux density for 
all spatial structure s. Details on the radio image can be found in 
iPetruket aLl(l2008l) (see their Sect. 3 and Fig. 3). In the fittings of 
the X-ray spectra the radio flux at 1 GHz is fixed at its best value 
(derived from the radio map), but we verified that allowing it to 
vary within its uncertainties (it is determined with a ~ 20% -50% 
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accuracy, depending on the spectral regions), the best-fit values 
and errors do not change significantly. The spectral index a is not 
fixed in our analysis, at variance with Rothenflug et al. (2004), 
who fix a = 0.6. 

We find relatively large residuals between data and spectral 
model in the 0.7-0.85 keV band. The same pro blem has been ob- 
served in the S UZAKU spectra of SN 1006 by Yamagu chi et al.l 
(I2008h and it has been interpreted as a result of the presence of 
K-shell O VII emission lines (transition series larger than Ky) 
not included in the spectral code. To take int o account the Ks, 
K^, an d O VII lines we follow the recipe bv I Yamaguchi et aTl 
(12008') and we add three narrow gaussians (a similar procedure 
has been adopted also by Acero et al. 2007). We fix the Nh pa - 
rameter to 7 x lO^** cm~^, in agreement with Du bner et al. I (l2002h . 

The determination of the chemical abundances is compli- 
cated by the combined presence of thermal and non-thermal con- 
tinuum. We focus only on the elements whose line complexes are 
visible in the spectra of the "thermal" regions (i. e. O, Ne, Mg, 
and Sfl see upper panel of Fig. |2|, but we find that, even in ther- 
mal regions, the abundances are not well constrained. Although 
the relative ratios of O, Ne, Mg, and Si abundances are quite 
well determined in each thermal region, their absolute best-fit 
values are entangled with the normalization of the VPSHOCK 
component. In particular, for a given set of best-fit abundances, 
we obtain similar values by multiplying all the abundances 
by the same factor k (even with k = 100) and by reducing the 
normalization by the same factor To avoid this degeneracy (that 
also hampers the determination of the emission measure), we 
consider two large regions: the first one is the sum of regions 28, 
29, 0, and 1 (region SE); the second one is the sum of regions 
13-16 (region NW). The larger statistics in these large regions al- 
low us to accurately determine the chemical abundances, whose 
best-fit values are shown in Table |2] Notice that in the SE region 
the chemical composition of the X-ray emitting plasma is quite 
different from that in the NW region, where we find a lower Si/O 
ratio. In both cases, however, we observe significant overabun- 
dances. Moreover, the large Si/O ratios obtained in SE and NW 
regions are typical of type la SNe. We therefore conclude that 
the thermal component is associated with shocked ejecta. We 
assume regions 8 - 22 to have the same abundances as region 
NW and the remaining 15 regions to have the same abundances 
as region SE. 



Table 2. Best-fit parameters. 



Element 


Region NW* 


Region SE* 


O/Oo 


5.0 ±0.1 


4.4 ± 0.3 


Ne/Neo 


4 7+0.1 

^•'-0.2 


1 5+°-2 




12 ±2 


15 ± 1 


Si/Si0 


29!? 


50 ±3 


kT (keV) 


0.41 ± 0.03 




EM'" (10" ) 


1.6 ±0.1 


1 53+°'2 


Tps (10'" s cm-3) 


1.6 ±0.2 




reduced (d. o. f.) 


1.85 (426) 


1.84 (234) 



** Union of regions 28, 29, 0, and 1 of Fig. [T] 
*** Emission measure per unit area. 



The results of our spectral analysis are summarized in Fig.|3j 



3.1.1. Non-thermal emission 

As shown in Fig. [3] we find that the azimuthal profile of the 
synchrotron cut-off frequency Vbreak pr esents a trend tha t is per - 
fectly consistent with that derived by iRothenflug et al.l (|2004|) . 
with the cut-off frequency increasing in the non-thermal limbs, 
although our break frequencies are systematically lower than 
their ones (by a factor of ~ 10). This discrepancy is due to 
the value of the radio spectral index a. In fact, we here derive 
a ~ 0.50 ± 0. 04 from our spec t ral fitt ings in the non-thermal 
regionfl while iRothenflug et al.l (|2004|) do not measure a, and 
assume a = 0.6. The lower values of the spectral index nat- 
urally implies lower cut-off frequencies. Notice that, as shown 
by the green points in Fig. [3] we obtain a ~ 0.5 also by using 
the "ESAS approach" for the analysis of the background, as de- 
scribed in Sect. |2l Our best-fit value a ~ 0.5 and our cut-off 
frequencies are in agreement with the recent results obtained by 
Allen etalJ(l2008h . 



3.1.2. Tliermal emission 

Figure [3] clearly shows that the emission measure per unit area, 
EM, of the thermal component significantly decreases in the 
bright non-thermal limbs (we obtain the same trend also analyz- 
ing the spectra with the ESAS method described in Sect.|2]i. In 
particular, at North-East (regions 19-23) and at South- West (re- 
gions 6-10), EM is a factor of S 4 lower than at North- West or 
South-East and in regions 9 and 19 it is consistent with zero. 
This is quite surprising because one would expect EM to be 
large in the non-thermal limbs, where the particle acceleration is 
supposed to be more efficient. In fact. lDecourchelle et alj (l2000h 
have shown that, if the main shock is modified, the post-shock 
density is larger, even behind the reflected shock (and this is our 
case, since we are observing the shocked ejecta). We investigate 
in detail the reason and the physical origin of this result in the 
next section. 

The temperature of the ejecta is not uniform and, in partic- 
ular, in the north-western part of the rim (regions 13-18 in Fig. 
[3j the temperatures are clearly lower. Unfortunately, in the non- 
thermal limbs the determination of temperature is affected by 
large errors, because of the low EM. 

The ionization time-scale Tps has been determined only 
in thermal regions and is quite uniform, with a value around 
10'" s cm"^ (in the non-thermal regions we have fixed it to 
its mean value in the adjacent thermal regions). We observe a 
large deviation from the average value only in region 13, where 
Tps = S.O^Qg X 10^ s cm"^. This is not surprising because re- 
gion 13 contains the bright oxygen knot studied in detail by 



IVink et alj (l2003h who found tnei ~ 2.35 + 0.07 x 10'^ s cm ^ 
Indeed lVink et al.l (l2003l) used the NEI model (while we use the 
VPSHOC K model) and their valu e should be compared with 
Tps/2 (see iBorkowski et"aLll200ll) . therefore our estimates are 
perfectly consistent wit h their on e s. No tice also that the relative 
abundances found by ,Vink et al.l ( 1200 3) in the oxygen knot of 
ejecta are quite si milar to those we derive for region SE, namely 
IVinketalJ(l2003h find 0:Ne:Mg:Si ~ 0.18:0.05:0.2:1 in the oxy- 
gen knot and we find 0:Ne:Mg:Si ~ 0.1:0.03:0.3:1 in region SE 
(see Table|2|l. This result further confirms that the thermal com- 
ponent is associated with the ejecta. 



^ We fixed the other abundances to the corresponding photospheric 
solar values, except for the S abundance that is tied to that of Si. 



^ In the thermal regions we therefore fix o- = 0.5, as shown in Fig. [3] 
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3.2. X-ray emitting ejecta in SN 1006 

Since the thermal emission is completely dominated by the 
shocked ejecta, it is interesting to study their spatial distribu- 
tion an d to obtain an image of their X-ray emission. Bamba et al. 
( |2008|) suggested that there are differences in the spatial distri- 
bution of thermal an d non-thermal emission in the north-eastern 
region of SN 1006. ICassam-Chenai' et al. I (120081) assumed that 
(in the south-eastern quarter of the shell) the X-ray emission in 
the "oxygen band" (0.5 - 0.8 keV) is completely associated with 
the ejecta. Indeed, even in this soft band, there is a non-negligible 
flux associated with the non-thermal component, specially in the 
non-thermal limbs (see Fig. |2]i. We use our results on the non- 
thermal component to estimate this contribution. 

The spectral properties of the SRCUT component are very 
robust and reliable and do not strongly depend on the model 
adopted for the thermal component, since the normalization of 
the synchrotron emission is derived from the radio data and the 
values of a and of the azimuthal profile of the roll-off frequency 
are in good agreement with those reported in literature. Using 
our result, we calculate that only in one third of the rim (regions 
28-2 and 18 - 22) the contribution of the thermal emission in 
the oxygen band is larger than the 80% of the total flux, while in 
fifteen regions (4-11 and 19 - 25) it is < 50%. This means that 
the 0.5 - 0.8 keV map of SN 1006 shown in the central panel of 
Fig.[T]cannot be used as a proxy of the ejecta. 

The procedure that we use to produce the image of the "pure" 
thermal emission in the oxygen band from our spectral results on 
the non-thermal emission is described in detail in Appendix A 
and is based on the assumption that the emission in the hard band 
(2 - 4.5 keV) is entirely non-thermal (for the validity of this as- 
sumption see Appendix A). Under this assumption, it is possible 
to use the spectral results to extrapolate the non-thermal emis- 
sion in the oxygen band from the hard band image (in the ap- 
pendix we show that the effects of the uncertainties in the best-fit 
parameters are negligible). Then, by subtracting the non-thermal 
contribution, we obtain the pure thermal image (i. e. the emis- 
sion of the shocked ejecta) in the 0.5-0.8 keV band. The map is 
shown in Fig. |4] and presents striking differences from the total 
image in the same band (central panel of Fig.lTJ, specially in the 
non-thermal limbs, where the emission of the ejecta is negligi- 
ble. 

In particular, in the spectral regions where the emission mea- 
sure of the thermal component is low (e. g. regions 8, 9, 19, 21, 
22, see Fig. [3]), the thermal image indicates that the filling factor, 
/, of the ejecta is at its minimum. Therefore, it is the low vol- 
ume filled by the X-ray emitting thermal plasma (in the spectral 
region) that causes the decreasing of the emission measure. In 
general, the comparison between Fig. [3] and Fig. |4] reveals that 
there is a clear relationship between / and EM: the variations 
of the emission measure of the thermal X-ray component can be 
associated with the different filling factors of the ejecta in the 
spectral regions. 

To verify that this conclusion is correct and that our map of 
the soft emission of the ejecta is reliable, we perform two tests: 
1. We consider 8 regions that, according to our thermal map, are 
expected to contain a larger filling factor of ejecta than regions 
19 - 23 and 6, 8, and 9. These new regions are indicated (in red) 
in the upper panel of Fig.[T] Figure |3] shows that the EM in these 
new regions (indicated by the red crosses) is significantly larger 
than that of the corresponding regions on the rim, as predicted 
by our thermal image. Notice that the extension along the line 
of sight in these new regions is only a factor ^ 10% larger than 
in regions 19 - 23 and 6, 8, and 9, therefore the much larger 



variations in EM are, indeed, due to larger filling factor of the 
ejecta. 

2. We exploit the large Chandra spatial resolution to extract 
Chandra ACIS-I spectra from the very narrow regions shown 
in Fig.|5] Region A covers a thin non-thermal filament, region B 
is immediately behind, where the emission in the 0.5 - 0.8 keV 
band is still large, but it is not associated with the ejecta, accord- 
ing to our thermal image, while region C covers a bright knot 
in our thermal image of the ejecta. The spectra extracted from 
these regions are shown in the lower panel of Fig. |5]and confirm 
the predictions of our thermal image: no O-line emission is vis- 
ible in the spectrum of region B, that is very similar to that of 
the non-thermal region A, while in region C the oxygen line is 
clearly visible. 

4. Discussion 

4. 1 . Properties of the ejecta 

Our spatially resolved spectral analysis shows that thermal X- 
ray emission can be completely associated with the ejecta and 
the large Si/O ratios obtained in SE and NW regions are typ- 
ical of type la SNe. Despite the overall spherical symmetry of 
the remnant, the distribution of the ejecta in the shell is not 
uniform. In particular, we detect inhomogeneities in the chemi- 
cal composition of the ejecta, with the south-eastern part of the 
rim characterized by lower Si abundance and larger Ne abun- 
dance than the north-western part. We also observe lower ejecta 
temperatures in the north-western part of the rim in agreement 
with the findings of A cero etaTI (l2007i) . The inter pretation of 
this re sult is not straightforward. As proposed by lAcero et al.l 
(I2007h . in the North- West the remnant is interactin g with a dense 
medium, originating a very bright H„ filament dWinkler et al.l 
2003), which is spatially correlated with the X-ray thermal emis- 
sion. The physical conditions are then different from that in the 
South-East, where the expansion of the shell occurs in a more 
uniform medium. It is possible that the different thermal condi- 
tions between NW and SE are due to different heating mecha- 
nisms: in particular we can suppose the ejecta to be heated by 
the reverse shock at SE and by a reflected shock, generated by 
the interaction of the main shock front with the density enhance- 
ment at NW. According to this scenario the ISM density should 
be lower behind the dense cloud, thus causing an almost free ex- 
pansion in that direction until the shock has reached the dense 
Ha cloud and generated the reflected shock. This may also ex- 
plain the differences in the abundances of the ejecta, since it is 
possible that the reflected and reverse shock originated at dif- 
ferent epochs (the interaction with the dense cloud is relatively 
recent) and have, therefore, interacted with different layers of 
the expanding ejecta. Nevertheless, according to this scenario, 
we would expect to find different ionization time-scales in the 
two regions, while we find that the ionization time-scales at NW 
are quite similar to those at SE (see Table|2]and Fig.O and this 
seems to indicate that the inhomogeneities in the abundances are 
intrinsic (i. e. are the result of an anisotropic explosion) and not 
due to differences between the reverse and the reflected shock. 

We can estimate the time elapsed after the shock impact, f, 
and the ejecta density n^j. We determine t from our best-fit val- 
ues of EM and Tps in the large NW and SE regions (Table |2]l. 
To this end we first derive the post-shock density of the ejecta 
from their emission measure per unit area. We first consider the 
NW and SE regions and, by assuming that the extension along 
the line of sight is L ~ 1.85 x 10'"^ cm (i. e. the length of the 
chord intercepted by the shell along the line of sight at the center 
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of these spectral regions), we derive the density iiej ~ 0.1/ ^[f 
cmT^ (where / is the filling factor of the ejecta) both at NW 
and SE. We perform the same estimate in each of the thermal 
regions (regions - 2, 13 - 17, and 27 - 29) to check for the 
presence of azimuthal variations of the density, but we do not 
find any evidence for it. We then derive the time elapsed after 
the shock impact in the NW and SE regions, by using the re- 
lationship: t ~ Tps Iriej ~ 2500 + 500 ^ff yr. Our value of t is 
reasonable, considering that / < 1 , because of the clumpy nature 
of the ejecta. 

4.2. The distance between contact discontinuity and biast 
wave sliock 

The image of the pure thermal emission of the ejecta, shown in 
Fig- HI allows us to trace the position of the contact discontinuity, 
CD, over the whole remnant. The distance of the CD from the 
blast wave position, BW, can be a powerful tracer for the pres- 
ence of efficient particle acceleration at the blast wave shock. 
iDecourchelle et all (l2000h have shown that by increasing the in- 
jection efficiency tj, the distance between BW and CD decreases. 
[Rothenflug et al. (2004) have shown that the non-thermal X-ray 
emission of SN 1006 is consistent with the polar cap scenario 
and this implies that 77 is larger in the non-thermal limbs, where, 
therefore, we expect the contact discontinuity to be closer to the 
shock front. COS measured the distance between BW and CD 
in the southeastern quadrant of the shell, finding that the ratio 
BW/CD decreases towards the non-thermal limbs, as expected, 
but reaching too low values (~ 1.0). Nevertheless, to derive the 
position of the contact discontinuity, they assumed the emission 
in the 0.5 - 0.8 keV band to be a tracer of the ejecta and we ver- 
ified that this assumption is valid only in a small portion of the 
shell. We therefore adopt the same diagnostics used by COS, but 
we here use our pure thermal image in the oxygen band (Fig.|4|i 
to trace the contact discontinuity over the shell. 

We consider the radial profiles of the surface brightness in 
our thermal image for 360 different angles (the azimuthal angle 
increases counterclockwise from the North direction) and, for 
each angle, we identify the CD as the position of the relative 
maximum in brightness at the largest distance from the centefl 
having a surface brightness larger than a given threshold. We fix 
the threshold at 10% of the flux of the brightest point in the pure 
thermal image. Since the pixel-size of the thermal map is ~ 50% 
of the EPIC point-spread function (FWHM), this procedure pro- 
vides a quite accurate estimate of the CD position. The associ- 
ated error is the addition in quadrature rule of the pixel size (8") 
and of the sigma used to smooth the image (one pixel). The posi- 
tion of the main shock front in the the rmal regions can be de rived 
from the deep Hg, image presented in I Winkler et al.l ( l2003h . cor- 
rected for the expansion of the remnant in the years elapsed from 
the optical observation to our observations (assuming an expan- 
sion speed of 4000 km/s). As explained by COS, it is difficult to 
identify the optical rim by using a quantitative criterion (e. g. a 
threshold in the radial profile), so we determine the location of 
the BW, inspecting by eye the optical image. Nevertheless, it is 
not possible to use the H„ image to trace the position of the blast 
wave shock in the non-thermal limbs, since the optical emission 
is quite noisy and does not present sharp features. Therefore, in 
these regions, we use the 2 - 4.5 ke V image (lower panel of Fig. 



The positions of the blast wave shock and of the contact disconti- 
nuity are measured with respect to the same center as that chosen by 
COS 



[TJ to trace the position of the blast wave shock, by adopting the 
same criterion as that used for the CD. 

Figure |6] shows the azimuthal profile of the ratio between 
BW and CD, rgw/CD, over the whole shell. In particular, to re- 
duce the noise level, we plot for each angle the moving aver- 
age of rsw/CD computed considering the ten adjacent points. By 
comparing our azimuthal profile of rgw/CD with that obtained in 
the southeastern quadrant by COS we notice that, although a re- 
markable agreement in the thermal regions is present, our rsw/CD 
does not drop to very low values near the non-thermal limbs (at 
odds with COS who finds rgw/CD ~ 1-0). This is because the 
pure thermal image in the oxygen band differs significantly from 
the total 0.5 - O.S keV image near the non-thermal limbs. The 
value of rsw/CD is quite uniform and does not present relevant az- 
imuthal variations except for a couple of regions (between ~ 40 
and ~ 60 degrees and between ~ 250 and ~ 260 degrees) where 
it reaches two maxima at ~ 1.15. We verify that in these two re- 
gions remarkable shock breakouts are present, as shown in Fig. 
Q We interpret these breakouts as the result of the propagation 
of the main shock in a locally rarefied medium. If the shock has 
recently encountered a low density inhomogeneity its speed is 
larger with respect to the adjacent regions, while the position of 
the reverse shock and of the contact discontinuity are not imme- 
diately affected. This explains the locally large values of rsw/CD- 
Elsewhere rsw/CD is limited in the range rsw/CD ~ 1.05 - 1.12. 



Table 3. Parameters of the 3-D MHD model 



ISM density 


0.05 cm-' 


Ambient magnetic field 


30 fiG 


Explosion energy 


1.3 X 10" erg 


Mass of the ejecta 


1.4 Me 


Shock speed after 1000 yr 


4750 km s"' 


Shell radius after 1000 yr 


8.65 pc 



We verify that this value is not consistent with that expected 
for a non-modified Sedov shock. To this end w e consider the 
3-D rn agnetohydrodynamic model developed by ' Orlando et al.l 
(I2007h to describe the expansion of a SNR through the magne- 
tized ISM, adopting a set of parameters appropriate for SN 1006 
(see Table O. The model includes no eventual magnetic field 
amplification, and no effects on shock dynamics due to back- 
reaction of accelerated cosmic rays. We explore different config- 
urations of the interstellar magnetic field by performing a set of 
numerical simulations u sing FLASH, an ad aptive mesh refine- 
ment multiphysics code dFrvxell et al.ll2006i) . In all the cases we 
find that Richtmyer-Meshkov (R-M) instability develops in the 
inner shell, as the forward and reverse shocks p rogress through 
the ISM and ejecta, respectively (e.g. iKane et al. 1999). After 
1000 yrs of evolution, the R-M instability determines a quite 
complex 3-D density and temperature structure of the contact 
discontinuity which extends a bit further with respect to the 
unidimensional case (where the instabilities are not described), 
leading to different values of the ratio rsw/CD- The results of 
our simulation at t=1000 yrs are shown in Fig. [8] the "solid" 
surface is a tracer of the ejecta material, while the density of 
the shocked ISM is shown in semi-transparent blue. In all the 
cases examined, we find that, according to our MHD model 
of a non-modified SNR shock, we expect rsw/CD > 1.16, at 
f - 1000 yrs. As shown in Fig. |6] we o bserve much lower 
values of rswiCD in SN 1006. As shown bv IDecourchelle et al.l 



M. Miceli etal.: SN 1006 



7 



( I2OOOI) and lEUison et alj (|2004|) . when efficient particle acceler- 
ation is present and a significant fraction of the shock energy is 
deposited in cosmic rays, the region between BW and CD be- 
comes narrower, thus making rgw/CD lower. Since we observe 
much lower values of rsw/CD in the whole rim with respect to 
the non-modified case, we conclude that in SN 1006 the shock 
is modified everywhere. 

If the ambient magnetic field were aligned in the SW-NE di- 
rection in the plane of the sky, with no component along the 
line of sight, we would expect lower values of rgw/CD in the 
non-thermal limbs, at odds with our findings. Nevertheless, if 
the aspect angle, (p, between the interstellar magnetic field and 
the line of sight is significantly different from a right angle (i. e. 
the interstellar magnetic field is not in the plane of the sky and 
((> ^ 10 degrees), we will not observe any significant variations 
of rsw/CD, since, in this case, we do not observe region where the 
acceleration process is more efficient edge-on. Therefore, our re- 
sults suggests that the aspect angle in SN 1006 is significantly 
different from 90°. A similar conclusion has been obtained by 
iPetruk et al.l (l2008l) by comparing observed and synthesized ra- 
dio maps of SN 1006. 

5. Summary and conclusions 

1 . The emission at the rim of the shell can be described by a 
mixture of thermal emission from plasma in non-equilibrium 
of ionization and non-thermal SRCUT emission. We com- 
pare two different approaches to X-ray spectra analysis, ob- 
taining perfectly consistent results with the two methods. 

2. We measure the radio-to-X-ray photon index of the non- 
thermal component, finding a ~ 0.5, and we confirm 
the azimuthal p r ofile o f the roll-off frequency found by 
iRothenflug et alj (|2004|) . revising the values of the frequen- 
cies in the light of the new value of a. 

3. Thermal emission can be associated with shocked ejecta. 
The temperature and the chemical composition of the ejecta 
are not uniform. In particular, we find lower temperatures 
and lower Si/O abundances at NW than at SB. 

4. Our pure thermal image in the 0.5 - 0.8 keV band shows 
that even in this soft band, the emission at the bright limbs is 
mainly associated with the non-thermal component. 

5. We derive the azimuthal profile of the ratio rsw/CD between 
the position of the blast wave shock and the contact discon- 
tinuity in the whole shell. rgw/CD is fairly uniform and its 
mean value is lower than that expected from a non-modified 
shock, even taking into account hydrodynamic instabilities 
that make rBw/CD lower than that derived from simple 1-D 
hydrodynamic models. We conclude that the shock is modi- 
fied everywhere and that the aspect angle is significantly dif- 
ferent from 90°. 
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Appendix A: Production of the pure thermal image 

We here describe the procedure we developed to produce the 
image of the thermal emission in the 0.5 - 0.8 keV band (we call 
this image TH). The scheme of the procedure consists in three 
steps: A) we assume the emission in the 2 - 4.5 keV band to be 
completely associated with the non-thermal component; B) we 
extrapolate the 2 - 4.5 keV image to produce an image of the 
non-thermal component in the 0.5 - 0.8 keV band, (we call this 
image MONTH) by using the results of our spatially resolved 
spectral analysis ; C) we subtract MONTH from the total image 
in the 0.5 - 0.8 keV band, thus obtaining the thermal image. In 
the following we describe in detail these three steps. 

A.1. Step A (assumption) 

We assume the thermal emission to be negligible in the 2 - 4.5 
keV band. Indeed the results of our spatially resolved spectral 
analysis show that this assumption is strictly valid in the non- 
thermal limbs (where the non-thermal flux is ^ 99% of the to- 
tal), while it does not hold in the thermal regions (e.g. regions 
28 - 2), where about 50% of the flux in the 2 - 4.5 keV is asso- 
ciated with the thermal component. Nevertheless, we can show 
that this problem does not generate significant effects. In fact, 
even if we overestimate the contribution of the synchrotron emis- 
sion in the thermal region (by a factor of a few) this contribution 
is still small in the 0.5 - 0.8 keV band: for example we verified 
(after step B) that it is ~ 10% and ~ 6% of the total flux in the 
oxygen band in region and region 16, respectively. This means 
that, when we produce TH, the flux in the thermal regions is 
not significantly underestimated, while that in the non-thermal 
region is correct. In conclusion we can say that our assumption 
holds where it is more necessary. 

A.2. Step B 

From the 2-4.5 keV image we can produce MOMTH by consid- 
ering a given spectral shape of the synchrotron emission (i. e. a 
given value of a and a given value of the roll-off frequency Vhreak) 
in each pixel. Since we found that in the whole rim the photon 
index is fairly uniform (see Fig.[3]l, we can assume that in each 
point of the image in the hard band the spectrum has a photon 
index a - 0.5. As for Vbreak, we found in Sect. |3.1| that it depends 
on the azimuthal angle 9 at the rim. We then interpolate the val- 
ues of Vbreak shown in Fig.[3]by using a Fourier series (up to the 
6th order) and we assume that the azimuthal trend of VbreakiQ) 
at the rim is given by this interpolating function at angular dis- 
tances from the center R > Ri„,{0), where Ri„,{9) - 11.7' for 
regions 13 - 17 and Rmt(0) - 13.5' elsewhere. As for the radial 
profile of Vbreak{R) for R < Rint{0), we follow iRothenflug et alj 
( 1200 4*) and we assume that: i) in the center {R < 7') the cut- 
off frequency is uniform and its value is equal to the minimum 
Vbreak found in the rim and ii) for T < R < Ri„,{0), the cut-off fre- 
quency increases exponentially with R, going from its minimum 
value (at the center) to the corresponding value at the rim. In 
this way we produce the map of Vbreak shown in Fig. lA.ll Since 
we know the spectral shape of the synchrotron emission in each 
pixel of the 2 - 4.5 keV image, we can produce the image of the 
non-thermal emission in the 0.5 - 0.8 keV band. This image is 
shown in the lower panel of Fig. lA.ll We investigate how the un- 
certainties in the best-fit parameters influence MOMTH. Figure 
|A.2| shows the ratio, p, of the non-thermal EPIC MOS count-rate 
in the 0.5 - 0.8 keV and 2 - 4.5 keV bands as a function of Vbreak, 
obtained with a - 0.5 (black curve), a - 0.45 (blue curve), and 



a — 0.55 (red curve). This ratio provides the factor that allows us 
to produce, for each pixel, MOMTH from the 2 - 4.5 keV map. 
Since in Figure |A2] the red curve, the blue curve, and the black 
curve are very similar, we can conclude that the uncertainties in 
the determination of a do not significantly alter our results. As 
for the effects of the uncertainties in Vbreak, we focus only on the 
non-thermal limbs, since in thermal limbs, as explained in Sect. 
lA.ll the non-thermal count-rate in the 0.5 - 0.8 keV band is neg- 
ligible. Figure IA.2I shows that in the non-thermal limbs (where 
Vhreak ^ 4x lO"' Hz and where the errors in Vhreak are quite small, 
see Fig. [3]), p does not significantly depend on Vhreak- Therefore, 
in these regions, the effects of the uncertainties in Vhreak are also 
negligible. 

A.3. Step C 

Once MOMTH has been produced we simply subtract it from the 
total image in the oxygen band, thus obtaining the pure thermal 
image in the 0.5 - 0.8 keV band shown in Fig.lD 
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Figure 2. Upper panel: pn spectrum (upper, in black) and 
MOSl,2 (lower, in red and green) spectra of region 29 of Fig. 
[T]with the corresponding best-fit model and residuals. The best- 
fit thermal H-non-thermal model described in Sect. 13. H is shown 
as continuous line, while individuals components are shown as 
dashed black lines for the pn spectrum only. Relevant emission 
lines are also indicated. Lower panel: same as upper panel for 
region 23. 



Figure 1. Upper panel: mosaiced count-rate images (MOS- 
equivalent counts per second per bin) of SN 1006 in the 0.8 - 2 
keV band. The bin size is 8" and the image is adaptively 
smoothed to a signal-to-noise ratio 10. The 30 regions selected 
for the spectral analysis of the rim (in white) and the 8 regions 
selected for the study of the ejecta emission measure (in red, see 
Sect. I3.2l l are superimposed. North is up and East is to the left. 
Central panel: same as upper panel in the 0.5 - 0.8 keV energy 
band. A contour level at 10% of the maximum is superimposed. 
Lower panel: same as upper panel in the 2-4.5 keV energy band. 



0.5-0.8 keV 
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Figure 4. Mosaiced pure thermal image (MOS -equivalent 
counts per second per bin) of SN 1006 in the 0.5-0.8 keV band. 
The contribution of the non-thermal emission in the soft energy 
band has been subtracted. The bin size is 8" and the image is 
adaptively smoothed to a signal-to-noise ratio 10. The same con- 
tour level shown in the central panel of Fig.[T]is overimposed for 
comparison. 



Figure 3. Best-fit parameters for the 30 regions shown in the up- 
per panel of Fig. [T](/n black) and for the 8 regions of the up- 
per panel of the same figure {in red). The model consists in a 
non-equilibrium of ionization isothermal component plus a syn- 
chrotron SRCUT model. The blue crosses indicate where the pa- 
rameters are fixed. The errors are at 90% confidence level. The 
chemical abundances are tied to those of region NW, for regions 
8 - 22, and to those of region SE, for regions 23-7 (see Table 
|2]l. For comparison, the best-fit values of a obtained by using 
the "ESAS approach" for the analysis of the background, as de- 
scribed in Sect.|2]are shown as green stars. 
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channel energy (keV) 

Figures. Upper panel: Close-up view of region 24 of Fig. [T] 
The XMM-Newton count-rate images in the hard band is shown, 
together with the Chandra ACIS-I image in the 0.5 - 0.8 keV 
band compared with our XMM-Newton thermal image in the 
same band (Fig.|4]i. The three boxes indicate the regions where 
we extracted the Chandra ACIS spectra. Lower panel: Chandra 
ACIS spectra extracted from the three regions shown in the up- 
per panel 




Angle 



Figure 6. Azimuthal profile of the ratio of the position of the 
blast wave shock to the contact discontinuity with correspond- 
ing error bars. In order to reduce the fluctuations, for each angle 
we plot the moving average of the ratio, computed considering 
the ten adjacent points. The position of the contact discontinu- 
ity is derived from the pure thermal image in the oxygen band 
shown in FigE) while the position of the blast wa ve shock is de- 
rived f rom the Ha image of SN 1006 presented in I Winkler et al.l 
( l2003 b (black curve, missing points mark locations where Ha 
filaments are not well defined) and from the hard X-ray image 
shown in the lower panel of Fig. [1] (in red). The angle is mea- 
sured counterclockwise from the North. 




Figure 7. Close-up view of the 2 - 4.5 keV band image of Fig. [T] 
showing two shock breakouts in the rim of SN 1006. We super- 
impose two lines indicating the azimuthal angles 40 and 60 {left 
panel) and 250 and 260 {right panel) 
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Figure A.l. Upper panel: map of the cut-ofF frequency derived 
according to the procedure described in Sect. IA.2I The 30 re- 
gions selected for the spectral analysis of the rim are indicated 
in black. Lower panel: image of the non-thermal emission in the 
0.5 - 0.8 keV band. 
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Figure A.2. Ratio of the EPIC MOS count-rate in the 0.5 - 0.8 

keV and 2 - 4.5 keV bands as a function of the break frequency, 
obtained by assuming an SRCUT model with spectral index a - 
0.45 (blue curve), a = 0.5 (black curve), and a = 0.55 (red 
curve) 



